The Nuclear Suppression Factor for Pt increasing towards the Kinematic Limit 
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The suppression of high transverse momentum (Pt) hadron production in ultrarelativistic heavy- 
ion (A-A) collisions as compared to the scaled expectation from proton-proton (p-p) collisions ex- 
pressed as the nuclear modification factor Raa is experimentally well established and can be traced 
back to interactions between the hard parton shower and the soft bulk matter. Physical intuition 
suggests that the medium modification should cease to be important when the hard scale of parton 
production is much larger than the typical momentum scale in the medium (e.g. the temperature T) 
and that consequently Raa{Pt) — > 1 for Pt 3> T. However, Raa is not a 'simple' observable, but 
rather results from many different mechanisms which influence what happens when Pt increases. 
In particular, Raa does not necessarily approch unity even if the hadron momentum is probed at 
the kinematic limit of the reaction. The aim of this work is to identify and discuss such mechanisms 
and to present different model expectations of what one would find if one could measure suppression 
out to the kinematic limit for hard hadron production and give predictions for the Pt dependence 
at both RHIC and LHC. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

Jet quenching, i.e. the energy loss of hard partons created 
in the first moments of an A-A collision due to interac- 
tions with the surrounding soft medium has long been 
regarded a promising tool to study properties of the soft 
bulk medium created along with the hard process 
The effect of the medium is apparent from a comparison 
of high Pt hadron observables measured in A-A collisions 
with the same observables in p-p collisions. The current 
range of such observables includes the suppression in sin- 
gle inclusive hard hadron spectra Raa 0> the suppres- 
sion of back-to-back correlations [1, 0] and single hadron 
suppression as a function of the emission angle with the 
reaction plane [10(. Recently also preliminary measure- 
ments of fully reconstructed jets have become available 
B 

Single hadron observables and back-to-back correlations 
above 6 GeV (where hadron production is dominated by 
hard processes) are well described in detailed model cal- 
culations using the concept of energy loss [T^ - [T3 |, i.e. 
under the assumption that the process can be described 
by a medium-induced shift of the leading parton energy 
by an amount AE where the probability of energy loss 
is governed by a distribution P(AE), followed by a frag- 
mentation process using vacuum fragmentation of a par- 
ton with the reduced energy. This can be cast into the 
form of a modified fragmentation function (MMFF). If 
the vacuum fragmentation function, i.e. the distribution 
of hadrons produced from a parton at fractional momen- 
tum z given a hadronization scale /1 2 is D(z, fi 2 ), then the 
MMFF given the medium induced energy loss probability 
P(AE) can be written as 



D{z,p?) = 



e D 
dAEP(AE) — 



l-AE/fi'^' 



1 - AE/E 



(1) 
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Beyond the leading parton approximation in which en- 
ergy loss and fragmentation factorize, one has to solve the 
full partonic shower evolution equations in the medium 
while assuming that the non-perturbative hadronization 
takes place outside the medium. At least for light sub- 
leading hadrons in a shower, factorizing hadronization 
from the medium-modified parton shower is a reason- 
able assumption at both RHIC and LHC kinematics. 
There are several calculations which utilize such medium- 
modified showers analytically [T^-|l3] ■ Recently, also 
Monte Carlo (MC) codes for in-medium shower evo- 
lution have become available [l8l - [23l ] which are based 
on MC shower simulations developed for hadronic col- 
lisions, such as PYTHIA H or HERWIG H- These 
have, unlike current analytical computations, full energy- 
momentum conservation enforced at each branching ver- 
tex. In these calculations the MMFF is obtained directly 
rather than from an expression like Eq. (TT|). 
So far, the different pictures for the parton-medium inter- 
action have been explored and compared with data over 
a rather finite kinematical window with Pt < 20 GeV. 
There is a widespread expectation that if the Pt range of 
the measurement could be extended, either at RHIC or at 
LHC, one would eventually observe the disappearance of 
the medium effect. The origin of this expectation is that 
the medium is able to modify the hard parton kinemat- 
ics at a typical scale set by its temperature T, whereas 
the parton dynamics takes place at a partonic hard scale 
Pt, and if pt S> T the hard kinematics should be esen- 
tially unchanged, which can be realized for large hadronic 
Pt- For example, in the case of the nuclear suppression 
factor, this expectation would imply that Raa{Pt) ap- 
proaches unity for Pt 3> T. It is the aim of this paper to 
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discuss the physics contained in the shape of Raa(Pt) 
and to present what current models, both based on the 
energy-loss concept and on the in-medium parton shower 
concept, predict for the shape of Raa at very large mo- 
menta at both RHIC and LHC. 



II. NUCLEAR SUPPRESSION IN THE 
ENERGY LOSS PICTURE 

A. Qualitative estimates 

A qualitative argument why Raa should increase with 
Pt can be made as follows: Parton spectra can be ap- 
proximated by a power law as dN/dpx — const. /p T 
where n rs 7. Assume that one can approximate the ef- 
fect of the medium by the mean value energy loss (AE) 
(for realistic energy loss models, this is not a good ap- 
proximation, as fluctuations around the mean turn out 
to be large). In this case, the energy loss shifts the 
spectrum. This can be described by the replacement 
Pt Pt + (AE) in the expression for the parton spec- 
trum. Raa{pt) can then be approximated by the ratio 
of the parton spectra before and after energy loss as 

(2) 

and it is easily seen that this expression approaches unity 
for px 3> (AE). However, it is not readily obvious under 
what conditions the limit is reached even if the medium 
properties are known. Parametrically, the medium tem- 
perature T governs both the medium density and the 
typical medium momentum scale, but the total energy 
loss represents the cumulative effect of the medium, i.e. 
a line integral of medium properties along the path of the 
partons, and furthermore the physics of medium-induced 
radiation is rather complicated, interference between dif- 
ferent radiation graphs play a significant role, and there- 
fore the mean energy loss is not simply ~ T . Thus, in 
realistic calculation the mean energy loss at RHIC con- 
ditions is (AE) « O(10) GeV even for T < 0.35 GeV 
[111 ], and hence it can be understood that current data 
are relatively far from the limit. 

There are five main points which may be raised against 
the approximation Eq. @: 

• the estimate holds for partons and does not take 
into account fragmentation: 

This, however, is not a crucial issue for the ques- 
tion at hand. The fragmentation function D(z, /i 2 ) 
is steeply falling with z and as a result, fragmen- 
tation processes at low z are preferred. However, 
for given momentum scale of the hadron spectrum, 
low z implies high parton momentum, and this 
is suppressed because the parton spectrum is also 
steeply falling with pt- As a result, there is some 
typical intermediate (z) (dependent on hadron and 



parton type) which relates hadron and parton mo- 
mentum, for quarks fragmenting into light hadrons 
at RHIC kinematics (z) « 0.5 — 0.7. This means 
that hadronic Raa is to first approximation sim- 
ply scaled by this factor as compared to partonic 
Raa- Fluctuations around the average tend to 
smear out structures in the partonic Raa through 
the hadronization process, but do not alter the 
shape of Raa(Pt) beyond that. Thus, qualitatively 
Eq. ([2|) holds also on the hadronic level. 

• the estimate does not distinguish between quarks 
and gluons: 

This is moderately important, as energy loss is ex- 
pected to be stronger for gluons by a factor 9/4 
(the ratio of the Casimir color factors). At low 
Pt, hadron production is driven by gluonic pro- 
cesses as gluons are copiously available in the low 
x region in parton distribution functions (PDFs) 
28 -32]. However, hadron production at higher Pt 
probes more and more in the high x region in the 
parton distributions, and eventually valence quark 
scattering dominates. The hadronic Raa should 
therefore show a rise from gluonic Raa to the larger 
value of quark Raa which corresponds to the tran- 
sition from gluon to quark-dominated hadropro- 
duction. As shown in [2y], this is likely to be the 
mechanism underlying the rising trend observed in 
Raa at RHIC. For asymptotically high energies, 
the mechanism is not relevant however, as this is 
always a quark dominated regime. 

• the estimate neglects fluctuations around the aver- 
age energy loss: 

In the presence of fluctuations, P(AE) can be writ- 
ten as the sum of three terms, corresponding to 
transmission witout energy loss, shift of the parton 
energy by finite energy loss or parton absorption as 

P(AE) = TS(AE) + S ■ P(AE) + A ■ 6{AE - E) (3) 

where P(AE) is a normalized probability distribu- 
tion such that T + S + A = 1. Inserting this form to 
average over Eq. (U) with the proper weights, one 
finds 

RAA^f + JdAE-S- P(AE) (l - " . (4) 

It follows that Raa obtained with this expression 
is always bounded by T from below (if a fraction 
of partons escapes unmodified, no amount of mod- 
ification to the rest will alter this) and by (1 — A) 
from above (if partons are absorbed independent of 
their energy, Raa will never approach unity). In 
many calculations, A is determined by the condi- 
tion that a parton is absorbed whenever its cal- 
culated energy loss exceeds its energy, i.e. A and 
S are dependent on the initial parton energy. In 
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FIG. 1: (Color online) The LO pQCD parton spectrum for 
\fs = 200 GeV compared with a power law fit to the region 
from 20 to 50 GeV. 



particular, in the ASW formalism [33j, the energy 
loss can be formally larger than the initial energy, 
since the formalism is derived for asymptotically 
high energies E — > oo and small energy of radi- 
ated gluons AE <C E, but is commonly applied to 
kinematic situations in which these conditions are 
not fulfilled. Raa at given px is then equal to the 
transmission term T plus a contribution which is 
proportional to the integral of P(AE) from zero 
up to the energy scale E max of the parton, seen 
through the filter of the steeply falling parton spec- 
trum. Thus, in the presence of fluctuations, Raa 
is dominated by fluctuations towards the low AE, 
and Raa grows with pt since E max grows linearly 
with pt- If P(AE) includes fluctuations up to a 
maximal energy loss AE max , then for pt 3> AE max 
the original argument made for constant energy loss 
applies and Raa aproaches unity. In practice this 
may not be observable - the energy loss probabil- 
ity for RHIC kinematics may be substantial up to 
scales of O(100) GeV [H|, i.e. of the order of the 
kinematic limit. 

• the pQCD parton spectrum is not a power law: 

While in a limited kinematic range, the pQCD par- 
ton spectrum is approximated well by a power law, 
at about t/s/4 the power law fit becomes a bad de- 
scription of the spectrum. This is shown in Fig. [T] 
Close to the kinematic limit at s/s/2, the parton 
spectrum falls very steeply. If one would attempt 
a local power law fit in this region, the region of 
validity for the fit would be small and n very large. 
One can readily see from Eq. ^ that even for 
Pt 3> (AE) Raa does not approach unity when at 
the same time n — > oo. In other words, close to the 
kinematic limit, even a small AE causes a massive 
suppression simply because there are no partons 
available higher up in the spectrum which could be 
shifted down. For this reason, close to the kine- 
matic limit Raa 1 can not be expected, rather 
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FIG. 2: (Color online) Raa(Pt) calculated with nuclear ini- 
tial state effects only, obtained from the NPDF set [3r3 | and 
the EKS98 [3l[ set of nuclear parton distributions calculated 
for the whole kinematic range at RHIC. 



(dependent on the details of modelling) , something 
like Raa — > T should be expected. 

However, note also that the validity of factoriza- 
tion into a hard process and a fragmentation func- 
tion has been assumed for hadron production up to 
the kinematic limit. This may not be true, Higher- 
Twist mechanisms like direct hadron production 
in the hard subprocess (in the context of heavy- 
ion collisions, see e.g. [271 ] may represent a differ- 
ent contribution which, due to color transparency, 
remains unaffected by the medium at all Pt and 
which may be significantly stronger than fragmen- 
tation close to the kinematic limit. This could be 
effectively absorbed into a modified coefficient T 
which however ceases to have a probabilistic inter- 
pretation. 

• the nuclear initial state effects have not been taken 
into account: 

The initial state nuclear effects, i.e. the difference 
in nucleon [H, [2i| and nuclear [30H32T ] PDFs, are 
often thought to be a small correction to the final 
state medium effects. Over a large kinematic range, 
that is quite true. However, as one approaches the 
kinematic limit and forces the distributions into the 
x — > 1 valence quark distributions, one probes the 
Fermi motion region in the nuclear parton distri- 
butions where the difference to nucleon PDFs is 
sizeable. 

In Fig. Raa(Pt) is shown for RHIC kinematics 
taking into account only the nuclear initial state 
effects with two different sets of nuclear PDFs, but 
no final state medium induced energy loss. It is 
readily apparent that over most of the kinematic 
range, Raa{Pt) ~ 1, but that there is a strong 
enhancement visible above 80 GeV. 
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B. Detailed calculation 

The detailed calculation of Raa in the energy loss mod- 
els presented here follows the Baier-Dokshitzer-Mueller- 
Peigne-Schiff (BDMPS) formalism for radiative energy 
loss [2J using quenching weights as introduced by Sal- 
gado and Wiedemann [3j], commonly referred to as the 
Armcsto-Salgado- Wiedemann (ASW) formalism. 
The probability density P(xo,yo) for finding a hard ver- 
tex at the transverse position ro = (xo,yo) and impact 
parameter b is given by the product of the nuclear profile 
functions as 



r A (r + b/2)r A (r -b/2) 
T AA {h) 



(5) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density pa(t,z) as T^(r) = 
J dzpA{Y,z) and TaaO 3 ) is the standard nuclear overlap 
function T A A{b) = d?sT A (s)T A (s - b). 
If the angle between outgoing parton and the reaction 
plane is 4>, the path of a given parton through the medium 
C(t), i.e. its trajectory £ as a function of proper medium 
evolution time r is determined in an eikonal approx- 
imation by its initial position ro and the angle ip as 
C(t) = (xq + rcos(0), yo + rsin(0)) where the parton is 
assumed to move with the speed of light c = 1 and the 
x-dircction is chosen to be in the reaction plane. The 
energy loss probability P(AE) pat h for this path can be 
obtained by evaluating the line integrals along the eikonal 
parton path 



cj c (r o ,0)=/ dCCq(0 and (qL) (r„, 4>) = / dC$(C) 
Jo Jo 

(6) 

with the relation 



g(C) =if-2-e 3/4 (C)(coshp-sinhpcosa) (7) 

assumed between the local transport coefficient 
(specifying the quenching power of the medium), the en- 
ergy density e and the local flow ra pid ity p with angle a 
between flow and parton trajectory [3J,[35|. The medium 
parameters e and p are obtained from a 2+1-d hydrody- 
namical simulation of bulk matter evolution [36[ , chosen 
to have the RHIC and the LHC medium described within 
the same framework. w c is the characteristic gluon fre- 
quency, setting the scale of the energy loss probability 
distribution, and (qL) is a measure of the path-length 
weighted by the local quenching power. The parame- 
ter K is seen as a tool to account for the uncertainty 
in the selection of the strong coupling a s and possible 
non-perturbative effects increasingthe quenching power 
of the medium (see discussion in [9]) and adjusted such 
that pionic Raa for central Au-Au collisions is described 
at one value of Pa- 



using the numerical results of [33| and the defini- 
tions above, the energy loss probability distribution 
given a parton trajectory can now be obtained as a 
function of the initial vertex and direction (ro,</>) as 
P(AE;uj c (r,(t)),R(r,(f)))p ath = P(AE) path for uj c and 
R = 2u>c/(qL). From the energy loss distribution given 
a single path, one can define the averaged energy loss 
probability distribution P(AE))t aa as 



(P(AE)) 



Taa ' 



1 

"2tt 



dx / dy Q P(x ,y a )P(AE) path . 



(8) 

The energy loss probability P(AE) pat h is derived in 
the limit of infinite parton energy [331 ] . however in the 
following the formalism is applied to finite kinemat- 
ics. In order to account for the finite energy E of 
the partons (P(AE))t aa is truncated at AE = E and 
6(AE - E) J|° dAE P( AE) is added to the truncated 
distribution to ensure proper normalization. The physics 
meaning of this correction is that all partons are consid- 
ered absorbed if their energy loss is formally larger than 
their initial energy. The momentum spectrum of hard 
partons is calculated in leading order perturbative Quan- 
tum Chromodynamics (LO pQCD) (explicit expressions 
are given in [9j and references therein). The medium- 
modified perturbative production of hadrons can then be 
computed from the expression 



d^ h+X = Y, d ^ f+X ®(P(^E))T A ; 



W^ h (z,p?) 

( 9 ) 

where da^T^^ +x is the partonic cross section for the 
inclusive production of a parton /, D^ h (z, p?) the vac- 
uum fragmentation function for the hadronization of a 
parton / into a hadron h with momentum fraction z and 
hadronization scale p and from this the nuclear modifi- 
cation factor Raa follows as 



Raa(pt,v) 



dN\ A /d P Tdy 
T AA {b)daPP/dp T dy' 



(10) 



III. NUCLEAR SUPPRESSION IN THE 
MEDIUM-MODIFIED SHOWER PICTURE 

A. Qualitative arguments 

In a medium-modified shower picture, the whole par- 
tonic in-medium evolution of a parton shower following 
a hard process is studied, leading to a modification of 
the fragmentation function (FF) which is more general 
than Eq. (fT|). In this framework, Raa — > 1 is realized 
if the MMFF becomes sufficiently similar to the vacuum 
FF. A qualitative argument why the MMFF should ap- 
proach the vacuum FF for Px 2> T can be made by 
considering for example the RAD (radiative energy loss) 
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scenario of the MC code YaJEM (Yet another Jet Energy- 
loss Model). This model is described in detail in Refs. 

mm. 

The parton shower developing from a highly virtual ini- 
tial hard parton in this model is described as a series 
of 1 — > 2 splittings a — > be where the virtuality scale de- 
creases in each splitting, i.e. Q a > Qb,Q c and the energy 
is shared among the daugther partons 6, c as Ei, — zE a 
and E c = (1 — z)E a . The splitting probabilities for a 
parton a in terms of Q a , E a are calculable in pQCD and 
the resulting shower is computed event by event in a MC 
framework. In the absence of a medium, the evolution of 
the shower is obtained using the PYSHOW routine [13] 
which is part of the PYTHIA package [II]. 
In the presence of a medium, the main assumption of Ya- 
JEM is that the parton kinematics or the splitting prob- 
ability is modified. In the RAD scenario, the relevant 
modification is a virtuality gain 
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FIG. 3: (Color online) Ratio of medium-modified over vacuum 
quark fragmentation function into charged hadrons obtained 
in YaJEM for a constant medium with 5 fm length and q = 2 
GeV 2 /fm. Shown are the results for different initial quark 



(11) energies E. 



of a parton during its lifetime through the interaction 
with the medium. In order to evaluate Eq. (II ip during 
the shower evolution, the momentum space variables of 
the shower evolution equations need to be linked with a 
spacetime position in the medium. This is done via the 
uncertainty relation for the average formation time as 



Eh 

Ql 



Eb 

Ql 



(12) 



and randomized splitting by splitting by sampling from 
the distribution 



P{n) = exp 



n 
(n) 



(13) 



The limit in which the medium modification is unimpor- 
tant is then given by Q 2 AQ 2 , i.e. if the influence of 
the medium on the parton virtuality is small compared 
with the virtuality itself, the evolution of the shower takes 
place as in vacuum. Note that there is always a kinemat- 
ical region in which the condition can never be fulfilled: 
The region z — > 1 in the fragmentation function repre- 
sents showers in which there has been essentially no split- 
ting. Since the initial virtuality determines the amount of 
branchings in the shower, this means one probes events 
in which the initial virtuality Qq is not (as in typcial 
events) of the order of the initial parton energy Eq, but 
rather Qq ~ mh where nih is a hadron mass. Since rrih 
is, at least for light hadrons, of the order of the medium 
temperature, Q 2 3> AQ 2 can not be fulfilled in the re- 
gion z fv 1 of the MMFF — here the medium effect is 
always visible. 

This is illustrated in Fig. [3J Here the ratio of 
medium-modified over vacuum fragmentation function 



jjq^h ^ ag obtained in YaJEM is shown for a con- 
stant medium for different initial partonic scales fi p = E. 
For a low initial scale of E = 20 GeV, one observes that 
the whole range between z — 0.2 and z = 1 is sup- 
pressed in the medium, whereas the region below z = 0.1 
shows enhancement due to the hadronization of the ad- 
ditional medium-induced radiation. For larger initial 
scales, the region of enhancement becomes confined to 
smaller and smaller z and the fragmentation function ra- 
tio approaches unity across a large range. However, in 
the region z &s 1 suppression due to the medium always 
persists as expected. 

As a consequence, one can expect Raa —> 1 f° r Pt ^> 
T (where AQ 2 is assumed to be parametrically 0(T 2 )) 
except near the kinematic limit Pt ~ where the 

region z fa 1 of the MMFF is probed and suppression is 
expected to persist. 



B. Detailed calculation 

The detailed computation of Raa within YaJEM is out- 
lined in (20| . It shares many steps with the computation 
within the energy loss picture as described above, in par- 
ticular the medium averaging procedure. 
The basic quantity to compute is the MMFF, given a 
path through the medium. Due to an approximate scal- 
ing law identified in (l9l ]. it is sufficient to compute the 
line integral 



(14) 



in the medium to obtain the full MMFF D MM (z, ii 2 ,() 
from a YaJEM simulation for a given eikonal path of 



6 



the shower-initiating parton, where /x^ is the partonic 
scale. The link between q and medium parameters is 
given as previously by Eq. Q, albeit with a different 
numerical value for K. The medium-averaged MMFF is 
then computed as 

{Dmm{z,^ p ))t aa = 

/>27T f>OC f>00 

— d(f> dx dy P(x ,y )D MM (z, ^,Q)- 

Z7r J J-oo J ~oo 

(15) 

From this, the medium-modified production of hadrons 
is obtained from 

d^ h+X =J2 d ^ f+X ® (D MM (z,^)) TAA (16) 
/ 

and finally Raa via Eq. (|10p . A crucial issue when com- 
puting Raa for a large momentum range is that YaJEM 
provides the MMFF for a given partonic scale whereas a 
factorized QCD expression like Eq. ([TBI) utilizes a frag- 
mentation function at givem hadronic scale. In previous 
publications [19 . [20| , the problem has been commented 
on, but not addressed, as the variation in momentum 
scale for current observables is not substantial. In this 
paper, the matching between partonic and hadronic scale 
is done as follows: 

For several partonic scales, (Dmm(z, I^))t A a lSl com- 
puted, and the exponent n of a power law fit to the par- 
ton spectrum at scale fi p is determined. The maximum 
of z u (Dmm(z, fJ.p))r A A corresponds to the most likely 
value z in the fragmentation process, and thus the par- 
tonic scale choice is best for a hadronic scale Pt — z[i v . 
The hadronic Raa is then computed by interpolation be- 
tween different optimal scale choices from runs with dif- 
ferent partonic scales. Finally, in the region Pt — > y/s/2, 
(Dmm(z, s/4))t aa is always the dominant contribution 
to hadron production. 

The matching procedure between hadronic and partonic 
scale choice also leads to a significant improvement in the 
description of Raa in the measured momentum range at 
RHIC as compared to previous results [H, [2(| • 

IV. RESULTS FOR RHIC 

The nuclear suppression factor for 200 AGeV central Au- 
Au collisions at RHIC, calculated both in the energy 
loss picture (represented by the ASW model) and the 
medium-modified shower picture (represented by the MC 
code YaJEM), is shown over the full kinematic range in 
Fig. H and compared with PHENIX data 7]. For the 
ASW calculation, the partonic Raa is also indicated sep- 
arately for quarks and gluons. 

Before discussing details of the plot, let us recapitulate 
the main differences between the energy loss picture as 
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- - ASW (quarks only) 
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FIG. 4: (Color online) The nuclear suppression factor Raa 
at RHIC across the full kinematic range in 10% central 200 
AGeV Au-Au collisions. Shown are PHENIX data Q , a calcu- 
lation in the energy loss picture (ASW) with quark and gluon 
result shown separately, and a calculation in the medium- 
modified parton shower picture (YaJEM). 

exemplified by the ASW model and the medium-modified 
shower picture as represented by YaJEM: 

• ASW is derived for infinite parton energy, hence 
P(AE) is independent of the initial parton energy 
and the only energy dependence arises from the pre- 
scription to assign contributions where AE > E 
into an absorption term A. In contrast, YaJEM is 
a finite energy framework where the MMFF explic- 
itly depends on the initial parton energy. In partic- 
ular, within ASW there is an energy-independent 
transmission probability T which bounds Raa from 
below. 

• In the energy loss picture, it is not specified what 
happens to the lost energy. In contrast, within Ya- 
JEM the energy lost from the leading shower par- 
tons is recovered explicitly in a low z enhancement 
of the MMFF. 

In Fig |U these differences are apparent as follows: In 
the lowest Pt region from 6 GeV and above, there is 
small rise of Raa with Pt observed in ASW which is not 
seen in YaJEM. As apparent from the comparison of the 
ASW result for pions to the result for quarks and gluons, 
the rise in this region in the ASW model is driven by 
the transition from a gluon-dominated to a quark dom- 
inated regime — the ASW hadronic result subsequently 
approaches the quark result for larger momenta. This 
transition is also present in YaJEM, however it is masked 
by the onset of the low Pt enhancement, which just starts 
to become significant below 6 GeV and corresponds to a 
decreasing trend of Raa with increasing Pt- As a result, 
the two opposing effects roughly cancel and the YaJEM 
result appears flatter than the ASW result between 6 and 
25 GeV. 

For higher Pt, there follows a region up to 50 GeV in 
which both the ASW and the YaJEM result descrease 
slightly. This can be traced back to the fact that the 
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pQCD spectrum is not a power law, and that local power 
law fits result in increasing n for higher px- The two 
curves run in parallel till sw 75 GeV, then the predictions 
of the two models are strikingly different. 
The ASW curve turns upward beyond Pr — 75 GeV. A 
comparison with Fig. [5] shows that this has nothing to 
do with the final state energy loss, but reflects the Fermi 
motion region in the nuclear PDFs. At the kinematic 
boundary, the curve finally turns over to reach the trans- 
mission probability T, as all shifts in the spectrum at the 
kinematic boundary result in substantial suppression and 
the only remaining contribution are unmodified partons. 
In contrast, the YaJEM result shows a strong suppression 
from 75 GeV to the kinematic limit. This corresponds to 
the region z — » I in the MMFF in which suppression 
was always observed in Fig. [3l regardless of the initial 
energy. This suppression is strong enough to mask the 
enhancement from the nuclear PDF. In contrast to the 
ASW model, YaJEM does not include an ^-independent 
transmission term, thus Raa becomes very small towards 
the kinematic limit. In this, the finite-energy nature of 
the suppression in YaJEM is apparent. Note that the 
YaJEM result cannot be computed all the way to the 
kinematic limit due to lack of statistics in the MC re- 
sults at z — > 1. 

It is also clear that there is no region throughout the 
whole kinematic range at RHIC in either model for which 
Raa — > 1 could be observed. 

V. RESULTS FOR LHC 

It is then natural to expect that Raa 1 could be real- 
ized by probing even higher momenta beyond the RHIC 
kinematic limit, e.g. by studying Raa at the LHC. How- 
ever, in going to collisions at larger y/s, not only the kine- 
matic limit is changed, but also the production of bulk 
matter is increased, i.e. higher y/s corresponds to a mod- 
ification of both hard probe and medium. There is how- 
ever reason to expect that eventually one will find a re- 
gion in which P T ^> T, AE max , yj AQ 2 and Raa ~^ 1 can 
be realized: The kinematic limit y/s/ 2 increases linearly 
with y/s. However, the medium density does not. There 
are different models which try to extrapolate how the ra- 
pidity density of produced matter increases with y/s. The 
Eskola-Kajantie-Ruuskanen-Tuominen (EKRT) model is 
among the models with the strongest predicted increase, 
and has the scaling ^ rj y/s°' 574 . Thus, the rapidity 
density of bulk matter increases significantly slower than 
the kinematic limit for increasing y/s. 
Although the medium lifetime may increase substantially 
with y/s as well, the more relevant scale is the transverse 
size of the medium, as high px partons move with the 
speed of light and will exit the medium once they reach 
its edge. However, the transverse size of the medium is 
approximately given by the overlap of the two nuclei, and 
hence to Oth approximation independent of y/s (beyond, 
there is of course the weak logarithmic growth of all total 





' 1 ' 


1 ' 




— ASW 






— ASW (gluons only) 






— ASW (quarks only) 






- YaJEM 














w 


If II 






1 \\ 


JJ 
f 


i , ! , 1 


\ \ 

1 



500 1000 1500 2000 2500 
P T [GeV] 



FIG. 5: (Color online) The nuclear suppression factor Raa at 
LHC across the full kinematic range in 10% central 5.5 ATeV 
Pb-Pb collisions. Shown are results in the energy loss picture 
(ASW) with the quark and gluon result shown separately, and 
a calculation in the medium-modified parton shower picture 
(YaJEM). 

cross sections with y/s). Thus, for asymptotically high 
energies, the integration limits for a line integral along 
the parton path through the medium will not grow ar- 
bitrary large, and the integrand, i.e. the density distri- 
bution, will be the main change. All these arguments 
indicate that Raa 1 can thus be realized for large y/s 
despite the increased production of bulk matter. The re- 
maining question if the LHC energy y/s — 5.5 ATeV is 
large enough. 

The result of the detailed calculation shown in Fig. [S] 
indicates that this is not the case. For this calculation, a 
hydrodynamical evolution based on an extrapolation of 
RHIC results using the EKRT saturation model [36| has 
been used to account for the increased medium density 
and lifetime. All other differences to the RHIC result 
are either plain kinematics, or can be traced back to the 
scale evolution of the MMFF. 

As far as the shape of Raa(Pt) is concerned, the LHC 
predictions for both ASW and YaJEM agree, however 
quantitatively they differ substantially. At the heart of 
this difference is that ASW is an infinite energy formalism 
in which the larger y/s of LHC as compared to RHIC is 
chiefly reflected in the harder slope of the parton spectra, 
but not directly in P(AE). In contrast, within YaJEM, 
in addition to the harder slope of the parton spectrum, 
there is an explicit scale evolution of the medium effect 
in the MMFF (see Fig. [3]) . Since both mechanisms tend 
to increase Raa, the combined effect of scale evolution 
and parton spectrum slope leads, all things considered, 
to less final state suppression in the YaJEM result. 
The shape of Raa(Pt) can be understood by the mech- 
anisms also observed in the RHIC case. The initial steep 
rise and subsequent flattening reflects the changing slope 
of the parton spectra. Note that the transition from 
gluon dominated to quark-dominated hadron production 
is not an issue over most of the LHC kinematical range. 
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The final enhancement above 2 TeV is again driven by 
the Fermi motion region in the nuclear PDFs. Unlike in 
the RHIC case, at LHC kinematics the suppression ob- 
tained from YaJEM for this region is not strong enough 
to mask the enhancement. Finally, close to the kinematic 
limit, a small Raa is obtained. 

These results indicate that there is no reason to expect 
that the limit Raa — > 1 can be observed even with LHC 
kinematics. However, the general trend for larger Raa 
observed in the transition from RHIC to LHC indicates 
that the limit could be reached for asymptotically high 
energies over a large kinematic range, however not close 
to the kinematic boundary. 

VI. DISCUSSION 

So far, the nuclear suppression factor Raa has been ob- 
served experimentally only in a very limited kinematical 
region. In this region, no strong Pt dependence has been 
observed. The main expectation of how Raa changes 
if observed over a larger kinematical range is that the 
suppression should eventually vanish and Raa approach 
unity. The results presented here show that this expres- 
sion is too simplistic. 

In particular, it is wrong to think of the shape of 
Raa{Pt) to be the result of any single cause. Instead, 
many effects, among them the slope change of the pQCD 
parton spectrum, the scale evolution of the medium mod- 
ification effect, the transition from gluon-dominatcd to 
quark-dominated hadron production and also the initial 
state nuclear effects all influence Raa(Pt) in a character- 
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istic way. Moreover, it is not sufficient to think of going 
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appears unlikely that the simple limit Raa ~* 1 for suffi- 
ciently high Pt can be reached even at LHC kinematics. 

These findings may largely be of little practical value 
due to the impossibility of reaching out to a substan- 
tial fraction of the kinematic limit experimentally. How- 
ever, theoretically they serve well to illustrate that even 
a hard probe observable like Raa is never 'simple' in 
the sense that it reflects directly tomographic properties 
of the medium, but rather that it is a convolution of 
many different effects which all need to be understood 
and discussed carefully. In particular, Raa cannot be 
interpreted as an observable reflecting properties of the 
medium causing a final state effect. The shape of the 
underlying parton spectrum or initial state effects are 
equally important to understand Raa- 
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